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Abstract 

In this paper, we design linear time algorithms to 
recognize and determine topological invariants such as 
the genus and homology groups in 3D. These properties 
can be used to identify patterns in 3D image recogni- 
tion. This has many and is expected to have more ap- 
plications in 3D medical image analysis. Our method 
is based on cubical images with direct adjacency, also 
called (6,26)-connectivity images in discrete geometry. 
According to the fact that there are only six types of lo- 
cal surface points in 3D and a discrete version of the 
well-known Gauss-Bonnett Theorem in differential ge- 
ometry, we first determine the genus of a closed ID- 
connected component (a closed digital surface). Then, 
we use Alexander duality to obtain the homology groups 
of a 3D object in 3D space. This idea can be extended to 
general simplicial decomposed manifolds or cell com- 
plexes in 3D. 



1. Introduction 

In recent years, there have been a great deal of new 
developments in applying topological tools to image 
analysis. In particular, computing topological invariants 
has been of great importance in understanding the shape 
of an arbitrary 2-dimensionall (2D) or 3-dimensional 
(3D) object [8 |. The most powerful invariant of these 
objects is the fundamental group [71- Unfortunately, 
fundamental groups are highly non-commutative and 
therefore difficult to work with. In fact, the general 
problem in determining whether two given groups are 
isomorphic is undecidable (meaning that there is no al- 
gorithm can solve the problem) lfl5l . For fundamen- 
tal groups of 3D objects, this problem is decidable but 
no practical algorithm has been found yet. As a result, 
homology groups have received the most attention be- 
cause their computations are more feasible and they still 
provide significant information about the shape of the 
object ||6l (9) 13). This leads to the motivating prob- 
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lem addressed in this paper: Given a 3D object in 3- 
dimensional Euclidean space R 3 , determine homology 
groups of the object in the most effective way by only 
analyzing the digitization of the object. 

The properties of homology groups have applica- 
tions in many areas of bioinformatics and image pro- 
cessing J9j. We particularly look at a set of points in 
3D digital space, and our purpose is to find homology 
groups of the data set. 

Many researchers have made significant contribu- 
tions in this area. For R 3 , based on simplicial decom- 
position, Dey and Guha have developed an algorithm 
for computing the homology group with generators in 
0(n 2 -g), where g is the maximum genus among all dis- 
connected boundary surfaces. This algorithm has been 
improved by Damiand et al 0]. They used the bound- 
ary information to simplify the process. However, the 
overall time complexity remained the same. 

In 2D, both R 2 and the cubical complex of linear 
algorithms (which is similar to that of digital spaces) are 
found to calculate Betti numbers, which are essentially 
the same as the genus |5| and |8|. 

Other algorithms for homology groups in cubical 
spaces are studies in 2D, these algorithms are either 
0(nlog 2 n) in lfl3l or 0(nlog 3 n) in (8). In general, the 
homology group can only be computed in 0(n 3 ) time 
for the cubical complex in (8). More about computa- 
tional homology is discussed in |9j . 

In this paper, we discuss the geometric and algebraic 
properties of manifolds in 3D digital spaces and the op- 
timal algorithms for calculating these properties. We 
consider digital manifolds as defined in [|3j|. More in- 
formation related to digital geometry and topology can 
be found in [11] and fl2l . 

In this paper, we introduce an optimal algorithm with 
time complexity 0{n) to compute genus and homology 
groups in 3D digital space, where n is the size of the 
input data. In Section 2, we review some properties of 
digital surfaces and manifolds [3|. Based on the clas- 
sical Gauss-Bonnet Theorem, we calculate the genus 
of a digital closed surface in 3D. Section 3 covers the 



necessary background in 3-manifold topology. Using 
Alexander duality, we relate homology groups of a 3D 
object to its 2-dimensional boundaries. In Section 4, we 
present our algorithm for homology groups. 

2 Gauss-Bonnet Theorem and Closed 
Digital Surfaces 

Cubical space with direct adjacency, or (6,26)- 
connectivity space, has the simplest topology in 3D dig- 
ital spaces. It is also believed to be sufficient for the 
topological property extraction of digital objects in 3D. 
Two points are said to be adjacent in (6,26)-connectivity 
space if the Euclidean distance of these two points is 1, 
i.e., direct adjacency. 

Let M be a closed (orientable) digital surface in the 
3D grid space in direct adjacency. We know that there 
are exactly 6-types of digital surface points [3]]||2]. 




Figure 1. Six types of digital surfaces 
points in 3D 

Assume that M; (M 3 , Af 4 , M 5 , M 6 ) is the set of 
digital points with i neighbors. We have the following 
result for a simply connected M J3) : 

|M 3 |=8 + |M 5 |+2|M 6 |. (1) 

M4 and Mq has two different types, respectively. 

The Gauss-Bonnet theorem states that if M is a 
closed manifold, then 

/ KadA — 2irx(M) (2) 

J M 

where dA is an element of area and Kg is the Gaussian 
curvature. 

Its discrete form is 

E { P is a point in m] K(p) = 2tt ■ (2 - 2g) (3) 
where g is the genus of M. 



Assume that Ki is the curvature of elements in Mi, 
i = 3,4,5,6. We have 

Lemma 2.1 (a) K 3 = it/2, (b) K 4 = 0, for both types 
of digital surface points, (c) K§ = —tt/2, and(d) Kq = 
—it, for both types of digital surface points. 

Proof. We can see that K4 is always since 
Kg = K\ x K2 and one of the principal curva- 
tures must be 0. We know that there exists a simply 
closed surface that contains only 8 points of M 3 
and several points within M4. Based on (1) and (3), 
8 • K 3 = 2tt • (2 - 0). Therefore, K 3 = tt/2. There 
is a simply closed surface that only contains M 3 , M4, 
and M 5 . So, K 3 = -K 5 . Same for K 6 = -2 ■ K 3 . So, 
K 6 = -2tt/2 = -7r. ■ 

Lemma 2.1 can also be calculated by the discrete 
Gaussian Curvature Theorem IfTBI . The curvature of 
the center point of the polyhedra is determined by 

/ K G dA = 2- ! r-Z i e i . (4) 

J M 

Gauss-Bonnet Theorem will hold based on (4). 

Therefore, we can obtain the same results as Lemma 
2.1. For example, in 3D digital space, the angle of one 
face is tt/2. So, K 5 = 2tt - 5 ■ tt/2 = -tt/2. (The 
proof of Lemma 2.1 is necessary if we are admitting 
the Gauss-Bonnet Theorem first and then obtaining the 
curvature for each surface point. We include a proof 
here since the reference [16] does not contain such a 
proof.) 

Given a closed 2D manifold, we can calculate the 
genus g by counting the number of points in M 3 , M5, 
and Mq. According to (3), we have 

Ef =3 /ir i .|M i |=27r-(2-2 fl ), 

tt/2 ■ \M 3 \ - tt/2 ■ |M 5 | -7T • |M 6 | = 2tt • (2 - 2g), 

|M 3 |-|M B |-2|M 6 |=47r-(2-2fl). 
Therefore, 

ff = 1 + (|Af 5 | + 2 • |M 6 | - |Af 3 |)/8. (5) 

Lemma 2.2 There is an algorithm that can calculate 
the genus of M in linear time. 

Proof. Scan through all points (vertices) in M and 
count the neighbors of each point. We can see that a 
point in M has 4 neighbors indicating that it is in M4 
as are M5 and Mq. Put points to each category of Mi. 
Then use formula (5) to calculate the genus g. ■ 



The two following examples show that the formula 
(5) is correct. The first example shown in Fig. 2, is 
the easiest case. In Fig. 2 (a), there are 8 points in A/3 
and no points in M5 or (To avoid the conflict be- 
tween a closest digital surface and a 3 -cell (3), we can 
insert some points on the surface but not at the cen- 
ter point.) According to (5), g = 0. Extend Fig. 2 (a) 
to a genus 1 surface as shown in Fig. 2 (b) where there 
are still 8 M3 points but 8 M5 points. Thus, Fig. 2 (b) 
satisfies equation (5). We can extend it to Fig. 2 (c), it 
has 16 M5 points and 8 M3 points. So g = 2. Using the 
same method, one can insert more handles. 

The second example came from the Alexander 
horned sphere. See Fig. 3. First we show a "U" shape 
base in Fig. 3(a). It is easy to see that there are 12 
A/3 points and 4 A/5 points. So g = according to 
Equation (5). Then we attach a handle to Fig. 3(a) 
shown in Fig. 3(b). We have added 4 A/3 points and 
12 M 5 points. g = l + (\M 5 \ + 2 • \M 6 \ - \M 3 \)/8 = 
1 + (4 + 12 - 12 - 4) = 1. Finally, we add another han- 
dle to the other side of the "U" shape in Fig. 3(a), the 
genus number increases by one since we still add 4 AT3 
points and 12 M5 points shown in Fig. 3(c). g = 2 for 
(c). For more complex cases like the Alexander horned 
sphere, one just needs to insert two smaller handles to 
an existing handle, so the genus will increase accord- 
ingly. 
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Figure 2. Simple examples of closed sur- 
faces with g = 0,1,2 



The above idea can be extended to simplicial cells 
(triangulation) or even general CW fc-cells. This is be- 
cause, for a closed discrete surface, we can calculate 
Gaussian curvature at each vertex point using formula 
(4). (The key is to calculate all angles separated by 1- 
cells at a vertex) Then use (3) to obtain the genus g. 
Since each line-cell (1-cell) is involved in exactly two 
2-cells, it is only associated with four angles. Therefore 
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Figure 3. An example came from Alexan- 
der horned sphere in digital space 



the total complexity will be 0{\E\) where E is the set 
of 1-cells (edges). Thus, 

Lemma 2.3 There is an algorithm that can calculate 
the genus of a closed simplicial surface inO(\E\) where 
E the set of 1-cells (edges). 

There are examples that \E\ is not linear to the num- 
ber of vertices \V\. 

3 Homology Groups of Manifolds in 3D 
Digital Space 

Homology groups are other invariants in topological 
classification. For a fc-manifold , Homology group Hi, 
i = 0, k indicates the number of holes in each i- 
skeleton of the manifold. Once we obtain the genus of 
a closed surface, we can then calculate the homology 
groups corresponding to its 3-dimensional manifold. 

Consider a compact 3-dimensional manifold in R 3 
whose boundary is represented by a surface. We show 
its homology groups can be expressed in terms of its 
boundary surface (Theorem |3.4t . This result follows 
from standard results in algebraic topology. Since it 
does not seem to be explicitly stated or proved in any 
standard reference, we include a self-contained proof 
here Q. This result follows from standard results in al- 
gebraic topology. It also appears in |6| in a somewhat 
different form. For the convenience of readers, we in- 
clude a short self-contained proof here. 

First, we recall some standard concepts and results in 
topology. Given a topological space Af, its homology 



groups, Hi(M), are certain measures of i-dimensional 
"holes" in M. For example, if M is the solid torus, 
its first homology group Hi(M) = Z, generated by 
its longitude which goes around the obvious hole. For a 
precise definition, see, e.g. Q. Let &j = rankHj(M, Z) 
be the ith Betti number of M. The Euler characteristic 
of M is defined by 

x(M) = j2(-iy% 

i>0 

If M is a 3-dimensional manifold, Hi(M) — for 
all i > 3 essentially because there are no i-dimensional 
holes. Therefore, x(M) = b a — bi + b 2 — &3- Further- 
more, if M is in R 3 , it must have nonempty boundary. 
This implies that 63 = 0. 

The following lemma is well known for 3 -manifolds. 
It holds, with the same proof, for any odd dimensional 
manifolds. 

Lemma 3.1 Let M be a compact orientable 3-manifold 
( which may or may not be in R 3 j. 

(a) IfM is closed (i.e. dM = %), then X {M) = 0. 

(b) In general, X {M) = §x(0M). 

Proof, (a) If dM = 0, the result follows from the 
Poincare duality which says that Hi(M) = H 3 ~ l (M), 
and the Universal Coefficient Theorem which says that 
the free part of H 3 ~ l (M) is isomorphic to the free part 
of i?3_j(M). Together, they imply that fcj = 63-^. 
Hence x(M) = b - b\ + b 2 - b 3 = 0. 

(b) In general, we consider the double of M 
denoted by DM, which is obtained by gluing two 
copies of M along dM via the identity map. By 
what we just proved, x(DM) = 0. On the other 
hand, the Euler characteristic satisfies a nice additive 
property: x(DM) = x( M ) + x(M) - xipM) = 
2 X (M) - x(dM). This implies X (M) = \x{dM). ■ 

Next, we recall the Alexander duality. 

Proposition 3.2 Let X C S n be a compact, lo- 
cally contractible subspace of S n where S n is the n- 
dimensional sphere. Then 

Hi(S n - X) = ir^-^X) for all i where H is the 
reduced homology. 

We remark that S n is the one point compactification 
of R n . Therefore, a submanifold in R n is automati- 
cally considered as a submanifold in S n in a natural 
way. Conversely, a submanifold M in S n is automati- 
cally a submanifold in S n unless M = S n . 

Before we prove Theorem 13 .41 we first prove a spe- 
cial case when dM is connected. 



Lemma 3.3 Let S be a closed connected surface in S 3 . 

(a) Its complement, S 3 — S, must have exactly two con- 
nected components. We denote them by M and M'. 

(b) ffi(M) ~ H X {M') Z^ s \ H 2 (M) = 
H 2 {M') = 0. 

Proof, (a) By the Alexander duality, H (S 3 - S) = 
H 2 (S) = Z. Therefore S 3 — S must have exactly two 
components. 

(b) Again, by the Alexander duality, 
H 2 (S 3 - S) H°(S) = 0, and this im- 

plies H 2 (M) = H 2 {M') = 0. We also have 
Hi(M) © Hi(M') = Hi(M U M') = H^S 3 - S) ^ 
H 1 (S) = Z hl ^ s \ In particular, there is no torsion in 
Hi (M) or H x (M'). By LemmaEQ 1 -rankffi (M) = 
1 - rankffi(M') = ± X (S) = |(2 - h(S)). This 
implies H X (M) ~ H X {M) ^ Z? bl ^ ■ 

Now we consider a general compact connected 3- 
manifold M in S 3 . Its boundary, dM, is a closed ori- 
entable 2-dimensional manifold possibly with several 
components. 

Theorem 3.4 Let M be a compact connected 3- 
manifold in S 3 . Then 

(a) H a {M) = Z. 

(b) Hi(M) £3 Z? bl(9M \ i.e. Hi(M) is torsion-free 
with rank being half of rank Hi (dM). 

(c) H 2 (M) = Z n ~ x where n is the number of compo- 
nents of dM. 

(d) H 3 (M) = unless M = S 3 . 

Proof. Statements (a) and (d) are obvious: (a) follows 
from connectedness of M, (d) is due to the fact that dM 
is non-empty. Next, we prove (c) and (d). 

Let Si , ■ ■ ■ , S n be the connected components of dM, 
By Lemma [3~3l each Si separates S 3 into two connected 
components, Mj and ML Since M is connected, it must 
be entirely contained in either Mj or ML Let M[ be 
the one containing M. It follows that S 3 — UiSi — 
M U UjMi. 

By the Alexander duality, H 2 (M U UjMj) = 
H 2 (S 3 - \JfSi) = H°(UiSi) S Z n -\ But H 2 (M{) = 
for each i (Lemma 13.3b . Therefore H 2 (M) = 
H 2 (MUU t Mi) S Z n ~\ 

Next, also by the Alexander duality, Hi(M U 
UiMi) = Hi(S 3 - UiSi) S H^UiSi) = H\dM) Si 
Z b ^ dM L But LHS = Hi(M) © ® i Hi(M i ) = 
Hi(M) © ®iZ^(Si) ^ Hi(M) © z^ dM L It 
follows that Hi(M) = Z? blidM \ ■ 



4 A Linear Algorithm of finding Homol- 
ogy Groups in 3D 

Based on the results we presented in Sections 2 and 
3, we now describe a linear algorithm for computing the 
homology group of 3D objects in 3D digital space. 

Assuming we only have a set of points in 3D. We 
can digitize this set into 3D digital spaces. There are 
two ways of doing so: (1) by treating each point as a 
cube-unit that is called the raster space, (2) by treating 
each point as a grid point. It is also called the point 
space. These two are dual spaces. Using the algorithm 
described in 0, we can determine whether the dig- 
itized set forms a 3D manifold in 3D space in direct 
adjacency for connectivity. The algorithm is in linear 
time. 

The more detailed considerations of recognition 
algorithms related to 3D manifolds can be found 
in Q] where Brimkov and Klette made extensive 
investigations in boundary tracking. The discussions of 
3D objects in raster space can be found in lfT4l . 

Algorithm 4.1 Let us assume that we have a connected 
M that is a 3D digital manifold in 3D. 

Step 1. Track the boundary of M, dM, which is a 
union of several closed surfaces. This algorithm 
only needs to scan though all the points in M to 
see if the point is linked to a point outside of M, 
That point will be on boundary. 

Step 2. Calculate the genus of each closed surface in 
dM using the method described in Section 2. We 
just need to count the number of neighbors on a 
surface, and put them in Mi, using the formula (5) 
to obtain g. 

Step 3. Using the Theorem 13.41 we can get Hg, Hi, 
H2, and H3. Ho is Z. For Hi, we need to get 
bi(dM) that is just the summation of the genus 
in all connected components in dM. (See |7| and 
0.) Hi is the number of components in dM. H3 
is trivial. 

Lemma 4.1 Algorithm 4.1 is a linear time algorithm. 

Proof. Step 1 uses linear time. We can first track 
all points in the object using breadth-first-search. We 
assume that the points in the object are marked as "1" 
and the others are marked as "0." Then, we test if a 
point in the object is adjacent to both "0" and "1" by 
using 26-adjacency for linking to "0." Such a point is 
called a boundary point. It takes linear time because 
the total number of adjacent points is only 26. Another 



algorithm is to test if each line cell on the boundary has 
exactly two parallel moves on the boundary [3|. This 
procedure only takes linear time for the total number of 
boundary points in most cases. 

Step 2 is also in linear time by Lemma 2.2. 

Step 3 is just a simple math calculation. For Ho, 
H2, and H3, they can be computed in constant time. 
For Hi, the counting process is at most linear. ■ 

Therefore, we can use linear time algorithms to cal- 
culate g and all homology groups for digital manifolds 
in 3D based on Lemma 2.2 and Lemma 4.1. 

Theorem 4.2 There is a linear time algorithm to cal- 
culate all homology groups for each type of manifolds 
in 3D. 

To some extent, researchers are also interested in 
space complexity that is regarded to running space 
needed beyond the input data. Our algorithms do not 
need to store the past information, the algorithms pre- 
sented in this paper are always Oilogn). Here, logn is 
the bits needed to represent a number n. 
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